Nonlocality of Kohn-Sham exchange-correlation fields in dielectrics 



OS 



in 

(N 



c3 



-a 
c 

o 
o 



> 
in 

(N 

in 
o 

Os 
-i— > 



i 

c 

o 
o 



X 



David Vanderbilt 

Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08855-0849 

(May 22, 1997) 

The theory of the macroscopic field appearing in the Kohn-Sham exchange-correlation potential for 
dielectric materials, as introduced by Gonze, Ghosez and Godby, is reexamined. It is shown that this 
Kohn-Sham field cannot be determined from a knowledge of the local state of the material (local 
crystal potential, electric field, and polarization) alone. Instead, it has an intrinsically nonlocal 
dependence on the global electrostatic configuration. For example, it vanishes in simple transverse 
configurations of a polarized dielectric, but not in longitudinal ones. 
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Ever since it became clear that electric polarization 
is indeed a well-defined bulk quantity in an insulating 
crystal E]-fl|, the status of this electric polarization in 
the Kohn-Sham (KS) density-functional theory (DFT) 
HU has become a topic of considerable interest |§,0-^2) . 
(Throughout this Letter, I refer to the exact version of 
DFT in which the true KS exchange-correlation (XC) 
functional is presumed known.) In Ref. |2j it was argued 
that the polarization of a crystalline insulator should be 
given exactly by DFT, on the basis that charge densi- 
ties are given exactly, and that any errors in polarization 
would show up as errors in charge densities at surfaces or 
interfaces. Gonze, Ghosez and Godby (GGG) [0J|] then 
pointed out that in order for the DFT polarization to be 
correct, the DFT XC potential would need to have a lin- 
ear spatial variation (i.e, an "XC electric field" £ xc would 
have to be present). They formulated a new version of 
DFT appropriate for crystalline insulators, in which the 
density n(r) and the electronic polarization P are shown 
to be uniquely related to the periodic part of the potential 
V(r) and the electric field £. This extended Hohenberg- 
Kohn (HK) || principle then allows the XC energy to 
be expressed as a functional of n(r) and P, instead of 
just n(r) alone. Recently, Martin and Ortiz ]12| have re- 
formulated and extended this analysis. While agreeing 
with many of the conclusions of GGG, they nevertheless 
appear to express some doubts about the GGG interpre- 
tation of the XC field, preferring instead to focus on the 
HK and KS descriptions for the change in polarization 
connected to a change in field. 

In this Letter, I present an analysis that clarifies the 
role of the XC field in the exact KS theory. I start by 
deducing the behavior of the XC potential for several 
simple configurations of a finite sample of spontaneously 
polarized dielectric material in vacuum. These examples 
illustrate misleading aspects of certain arguments given 
by the previous Refs. |2|,||,[l2|] . Briefly, it is now under- 
stood |?],||,[l2| that the local periodic charge density in 
some small region of the sample can be generated by any 
of a continuous family of KS potentials labeled by the 
choice of effective field £ c g = £ + £ xc or, equivalently, by 



the choice of electric polarization P c g , in the same small 
region of the corresponding fictitious KS system. Based 
on Refs. Q and S], one might assume that the correct 
choice would be the one that makes the polarization cor- 
rect, P c ff = P; while in a naive approach one would make 
the choice £ e g = £, at least for the case £ = 0. Here, I 
show that neither of the above choices is generally correct. 
Instead, the correct choice is inherently nonlocal, and de- 
pends upon the electrostatic configuration of the entire 
system. For example, for configurations in which P(r) is 
essentially longitudinal, the correct choice is P e g = P; 
but if essentially transverse, then £ xc = 0; and for more 
complicated geometries, neither simple choice is correct. 
This ultra-nonlocality of the XC potential appears to be 
an inherent complicating feature of the exact DFT the- 
ory. 

I begin by establishing some notation and reviewing 
some basic results of Refs. J7|8p^]. Consider a periodic 
insulating crystal with fixed lattice vectors specifying the 
unit cell, and an external electron potential consisting 
of a periodic and a linear part, V(r) = V(r) — e£ ■ r. 
A tilde, as on V , will be used to indicate a quantity 
having the periodicity of the unit cell, and £ is a uniform 
electric field. As long as £ is not too large, one can 
with very good precision identify a physical (although, 
strictly, metastable) state of the system, having periodic 
density, that is connected to the £ = ground state 
by slow adiabatic switching of £ 12-14). Letting n be 
this periodic density, we can then search for the non- 
interacting KS system for which the effective potential 
has the same linear part (same field £) but periodic part 
V^ s . These relations are those of the conventional KS 
theory applied naively to the periodic system, and can 
be summarized as 



{V,£} ^{n,£} 
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(here T and 'NT indicate 'interacting' and 'noninteract- 
ing' respectively). Alternatively, one can identify the 
electronic polarization P of the true interacting system 
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FIG. 1. Sample geometries of a spontaneously polarized 
material. Details are given in text, (a) Cubic sample with 
spontaneous polarization along +£, partially reduced by de- 
polarization fields arising from surface charges, (b) Slab ge- 
ometry in which surface charges P tot ■ h are precisely can- 
celled by appropriately chosen external planar charges ^fao. 
(c) Geometry composed of four domains in such a way that 
V-P = 0. 



and search for the non-interacting system that correctly 
reproduces both n and P; that is, 



{V,S} — {H,P} «— >{^,£efF} 
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This is the approach of GGG Unless dictated by 

symmetry as for a centrosymmetric crystal in zero field, 
there is no reason to expect £ e g = £, any more than 
we expect some particular Fourier component of V e g to 
match that of V. Thus, the "exchange-correlation field" 
defined as £ X c = £eS — £ is generally non-zero. (It is un- 
derstood that £ X c is not a true electric field, since it acts 
only on the electrons.) In either case, Eq. (1) or (2), it 
has to be supposed that the non-interacting system is 
also an insulator p5[ , and that the field £ or £ c g acting 
on this non-interacting insulator is again small enough so 
that a metastable state is well defined [Q . 

Either procedure, Eq. (1) or (2), is perfectly sensible 
in the abstract, but we now have to investigate whether 
and how it might apply to the case of a more realistic 
non-uniform configuration of a dielectric material. Con- 
sider, for example, the three geometries sketched in Fig. 
1; we consider first the electrostatic configuration of the 
physical (interacting) system for each case. Figure 1(a) 
shows a cubic sample of a spontaneously polarized ma- 
terial. For definiteness, let us take this material to be 
BaTiOa in the cubic perovskite structure with the atomic 
coordinates frozen as follows: the unit cell is ideal cubic, 
the Ba nuclei lie at the cube corners, the O nuclei lie on 
the cube face centers, and the Ti nuclei are displaced by 
a constant distance 5xi = 0.05 A along x from the cube 
centers. This material retains a gap of several eV and has 



a spontaneous (zero-field) polarization P 
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with both components P; on (nuclear plus core) and p(°J 
(valence electronic) lying along +x [[l6| . The geometry of 
Fig. 1(a) is such that the surface discontinuities in P to t(r) 
give rise to macroscopic surface charges on the left and 



right faces (jT7| , generating an electric field that partially 
depolarizes the sample. Thus, the electronic polarization 
at the point A in Fig. 1(a) will be somewhat reduced from 
p(°) . In the configuration of Fig. 1(b), this is avoided by 
placing infinitely thin sheets of additional external charge 
=F<7 = ±Ptol on the left and right faces respectively, pre- 
cisely canceling the depolarization field Thus, the 
macroscopic electric field vanishes at point B, and the 
electronic polarization at this point is just P = P^x. 
Finally, Fig. 1(c) shows a different configuration in which 
the internal fields also vanish. This time the sample is 
comprised of four domains in which the displacements 
of the Ti nuclei are along e(r) with e = +x, —y, — x, 
and +y in the top, right, bottom, and left domains, re- 
spectively. Clearly the solution is £ (r) = since then 
P(r) = P(°>g(r) and thus V • P tot = V • (P ion + P) = 0, 
which is consistent with £ = 0. Thus, the local condi- 
tions at point B of Fig. 1(b) and point C of Fig. 1(c) are 
identical: £ = and P = P^x. 

Of course it must be assumed that the samples in Fig. 1 
are sufficiently large that macroscopic fields can be de- 
fined. So, when we speak of "point A," we really refer to 
a region, large compared to atomic dimensions but small 
compared to sample dimensions, in which a periodic V 
and field £ can be identified. But note, also, that for 
a macroscopic sample having the configuration of Fig. 
1(a), the electrostatic potential difference between the 
left and right faces may greatly exceed the band gap, so 
that in principle the ground state would become metallic. 
However, in the context of dielectric theory one is again 
much more interested in the metastable ]T^ ] insulating 
state obtained by starting from a configuration without 
macroscopic electric fields, such as that of Fig. 1(b), and 
then adiabatically restoring the fields. Throughout this 
Letter it will always be assumed that both the physi- 
cal and the fictitious KS systems are in such metastable 
states Q. 

Let us now deduce what must be the behavior of V xc (r) 
for each of the configurations of Fig. 1. The dot labeled 
'A' in Fig. 2 represents the values of the physical electric 
field £ and polarization P of the interacting system of 
Fig. 1(a) at point A. These determine the density n at 
point A, which must be reproduced by the fictitious KS 
system at point A. The dashed curve indicates the locus 
of values (£ e g, P e ff) of the KS system that are consistent 
with each other and with this given n. The choice of 
Eq. (1) corresponds to point A' (insisting that £ c g = £), 
while that of Eq. (2) corresponds to point A" (insisting 
that P c ff = P) ■ Any point on the dashed curve generates 
the correct periodic density at A, and so is a candidate 
for the state of the KS system at A. 

Now comes the crucial point of the argument. In order 
to decide which point on this curve should be selected, 
it is necessary to inspect the configuration of the system 
as a whole. This is illustrated by the geometries of Fig. 



2 




8 



FIG. 2. A, B, and C indicate the local state of the physical 
(interacting) system at corresponding points of Fig. 1; solid 
line is the dielectric relation for the physical material. A' 
and C indicate states of the Kohn-Sham system having the 
correct density and field, Eq. (1); A" and B" indicate KS 
states having the correct density and polarization, Eq. (2). 
Dashed and dotted curves indicate the KS dielectric relations 
defined under the constraint that the periodic density match 
that of the physical system at A, or B and C, respectively. 

1(b) and (c), for each of which the correct choice is easily 
deduced. As remarked above, in both cases the physical 
system is locally the same (point labeled 'B,C'). For the 
geometry of Fig. 1(b), the translational symmetry |ll| 
along y and z insures that both P and P e ff lie along 
x. But in either the physical (interacting) or the KS 
systems, the value of P is necessarily related || to the 
presence of a macroscopic electronic surface charge P ■ n; 
and since this electronic charge must be given correctly 
in the KS theory, it is safe to conclude that P c ff = P. 
Thus, for every point in the interior of Fig. 1(b), the 
state of the KS system is given by B" in Fig. 2. On 
the other hand, in Fig. 1(c) both the physical and KS 
electric fields must vanish everywhere. Heuristically, this 
is suggested by the invariance of V e g(r) under fourfold 
rotation. More precisely, we can check that a consistent 
solution is obtained when, for every point in the interior 
of Fig. 1(c), the state of the KS system is given by C in 
Fig. 2. In this case the polarization is given incorrectly by 
the KS theory at every interior point. Nevertheless, since 
V • P e ff = everywhere, such an error is still consistent 
with the KS system having exactly the correct n(r). 

The example of Fig. 1(b) and (c) demonstrates that a 
knowledge of the local state of the physical interacting 
system (i.e., a knowledge of the material, and of £ or 
equivalently of P) is inherently insufficient to determine 
the corresponding local state of the KS system. Thus, a 
field £ X c exists at point B, but not at point C, in spite of 
the fact that the local states of the physical systems at 
B and C are identical. 

In general, the correct correspondence to the KS sys- 
tem must be determined by an analysis of the electro- 



static configuration as a whole, and need not reduce to 
either simple choice (£ Q g = £ or P c ff = P). For example, 
for point A of Fig. 1(a), the correct choice of KS sys- 
tem might correspond to any of the points on the dashed 
curve of Fig. 2. To identify the correct point, an analysis 
of the following type would need to be carried through. 
Assuming that the physical £(r) and P(r) are known: 
(i) At each point in space, obtain n from Eq. (1) or (2), 
and also obtain the dielectric relation P c ff = Pa(£ c s) 
for a non-interacting system constrained to have fixed 
n as £ e g is varied (e.g., the dashed or dotted curve in 
Fig. 2). (ii) Taking first £ xc — (i.e., £ G ff = £) ev- 
erywhere, compute the macroscopic charge density error 
Sn(r) = V- [P c ff(r) — P(r)]. (iii) Find a macroscopic field 
£xc(r) such that, when acting on the system of nonlinear 
dielectric materials specified by dielectric responses Vn, 
it induces precisely a cancelling charge density — Sn(r). 

The above is rarely a simple procedure in practice, but 
it illustrates the conditions under which a field £ xc is 
needed. In Fig. 1(a) and (b), one finds Sn ^ 0, so the 
XC field must be present; while for Fig. 1(c), Sn = 
and so £ xc = 0. In fact, it is evident that for a purely 
longitudinal configuration such as Fig. 1(b), Eq. (2) ap- 
plies and the KS field £ e s fails to match £; while for a 
purely transverse configuration as in Fig. 1(c), Eq. (1) ap- 
plies and instead the KS polarization P c fr fails to match 
P. It is easy to see that infrared- active longitudinal and 
transverse phonons in a polar insulator behave in a way 
analogous to Figs. 1(b) and (c), respectively. 

A few comments may be in order, (i) According to 
the theory of Refs. PHI], the polarization P is actually 
only well-defined modulo eR/$l, where R is a lattice vec- 
tor and SI is the cell volume. For the puposes of the 
arguments given here, we assume that the polarization 
differences [e.g., between domains in Fig. 1(c)] are suf- 
ficiently small that there is no difficulty in choosing the 
correct branch of P. (ii) When constructing the KS po- 
tential for the configurations of Fig. 1, attention would 
also have to be given to the microscopic details of the 
choice of V^ c (r) at surfaces and interfaces. However, local 
modifications to V xc (r) at a surface or interface will not 
affect the macroscopic electronic charge appearing there 
so that the arguments about macroscopic fields will 
not be affected, (iii) In all of the arguments given here, 
the lattice coordinates have remained clamped. A realis- 
tic theory of dielectric materials would of course have to 
take into account the lattice relaxation effects. 

In summary, it has been shown that the exact KS po- 
tential has an ultra-nonlocal dependence upon the elec- 
tronic charge density, and this dependence does not take 
such a form that it can be automatically simplified by 
reference to the local electronic polarization. In some re- 
spects, the present conclusions are disappointing. One 
potentially appealing motivation for introducing [j^] 
the expanded HK correspondence {V, £} «-» {n, P}, ex- 
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pressed explicitly in Ref. JlJ] , is a hope that the general- 
ized XC functional of density and polarization might be 
more localized than the conventional functional of den- 
sity alone. The present conclusions unfortunately do not 
support this view. In fact, they reinforce previous indi- 
cations that the true KS functional V^ c [n] is, in general, 
an extraordinarily non-local and irregular functional of 
the density [|ll) . As for the peculiar properties of V xc [n] 
demonstrated here, there seems to be little hope of cap- 
turing these in practical approximations, even semilocal 
ones such as the weighted- or average-density approxi- 
mations |2(J] . And indeed, it is not obvious that it would 
be desirable to do so. Of course, the approximate local 
and generalized-gradient functionals currently in use do 
not have these peculiarities. 

This work was supported by NSF Grant DMR-96- 
13648. I wish to thank R.M. Martin for a conversation 
that stimulated the present work. 
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